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ASTRONOMY 

AND 
ASTROPHYSICS 



Abstract. We present results from observations obtained 
with ISOPHOT, on board the ISO satellite^, of a represen- 
tative sample of seventeen CSS/GPS radio galaxies and of 
a control sample of sixteen extended radio galaxies span- 
ning similar ranges in redshift (0.2 < z < 0.8) and radio 
luminosity (P2.7GHZ > 10 26 W/Hz). The observations have 
been performed at A = 60, 90, 174 and 200 fxm. The orig- 
inal purpose of these observations was to check whether 
CSS/GPS sources are associated with very gas rich galax- 
ies, as required by the scenario in which the growth of the 
radio source is inhibited by the dense medium of the host 
galaxy. 

Unfortunately the resulting performance of ISOPHOT 
was worse than expected. As a consequence, the detection 
limit at 60/im is similar to that obtained previously with 
IRAS but better than that at 90/im. 

Seven of the CSS/GPS sources have detections > 3cr at 
one or more wavelengths, one of which is detected at > 5a. 
For the comparison sample five objects have detections 
> 3cr one of which is at > 5a. 

By co-adding the data we have obtained average flux 
densities at the four wavelengths. 



Send offprint requests to: F. Pozzi 
Correspondence to: cfanti@astbol.bo.cnr.it 

1 Based on observations with ISO, an ESA project with in- 
struments funded by ESA Member States (especially the PI 
countries: France, Germany, the Netherlands and the United 
Kingdom) with the participation of ISAS and NASA. 



We found no evidence that the FIR luminosities of 
the CSS/GPS sources are significantly different from those 
of the extended objects and therefore there is not any 
support for CSS/GPS sources being objects "frustrated" 
by an abnormally dense ambient medium. 

The two samples were then combined, providing FIR 
information on a new sample of radio galaxies at in- 
termediate redshifts. We compare this information with 
what previously known from IRAS and discuss the aver- 
age properties of radio galaxies in the redshift range 0.2 - 
0.8. The FIR emission cannot be accounted for by extrapo- 
lation of the synchrotron radio spectrum and we attribute 
it to thermal dust emission. The average FIR luminosity 
is > 6 x lO 11 ^©. Over the observed frequency range the 
infrared spectrum can be described by a power law with 
spectral index a ~ 1.0 ± 0.2. Assuming the emission to 
be due to dust, a range of temperatures is required, from 
> 80 K to pa 25 K. The dust masses required to explain 
the FIR emission range from 5 x 10 5 M Q for the hotter 
component up to 2 x 10 s Mq for the colder one. 

We present also observations on four nearby (z < 0.1) 
GPS radio galaxies, two of which are detected at all four 
wavelengths. 

Key words: Galaxies: active; Galaxies: ISM; Infrared: 
galaxies; Infrared: ISM; continuum; Radio continuum: 
galaxies 
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FIR properties of CSS/GPS radio galaxies 



1. Introduction 

Compact Steep Spectrum (CSS) and GHz Peaked Spec- 
trum (GPS) radio sources are powerful extragalactic radio 
sources with radio emission confined well within their host 
galaxy/quasar (<^ 15kpc)^J They are as powerful as the 
classical FPJI radio sources but are of much smaller phys- 
ical size and yet have normal/steep radio spectra at GHz 
frequencies. A discuss ion of the p roperties of CSS can be 
found in Fanti et al. ( 1990| ; 1995). A comprehensive and 
updated review on CSS and GPS is presented by O'Dea 

Some 15 to 25 % of sources in a flux density lim- 
ited sample, depending on the selection frequency, be- 
long to this class. Many of them have radio spectra which 
show a flattening or a more marked turnover at fre- 
quencies between 50 MHz and a few GHz; those sources 
with a turnover at about 1 GHz are the GPS sources. 
The turnover is usually attributed to synchrotron self- 
absorption, but also free-free absorption (Bicknell et al. 



1997 ) and induced Compton scattering (Kuncic et al. 



1998) have been considered. 



Given that CSS and GPS sources are physically small, 
i.e. of sub-galactic dimensions, it has been suggested that 
they are either (1) young objects (youth scenario, Phillips 
& Mutel 1982| ), which have yet to develop extended radio 
lobes or (2) sources where the radio emitting plasma is 
trapped by an unusually dense interstellar medium (frus- 
tration scenario, van Breugel et al. 1984). 

In the latter hypothesis, CSS/GPS radio sources would 
be as old as those of larger size but their jets would spend 
their whole lifetime trying to escape, without success, out 
of the interstellar medium. This scenario requires a rather 
dense interstellar medium (average density 7Jism ^> 1 
cm~ 3 and total mass within 1 kpc A/ism > 1O 8 M ; see 



De Young 1993 



1995) which is in the path 



Fanti et al. 

of the radio source. Even if the host galaxy contains a 
lot of gas, if this is distributed in a disk perpendicular to 
the radio source axis, it will have little effect on the radio 
source propagation, therefore any "frustrating" gas has to 
be distributed over a large volume. Attempts to probe the 
different phases of this interstellar gas (e.g. optical obser- 
vations of the Narrow Lines, polarization studies, X-ray 
emission) have given, in our opinion, little support to the 
frustration scenario. 

Recently, proper motion measurements of the hot spots 
in a few GPS sources by Owsianik et al. (1998) and 



Owsianik & Conway (1998) have shown separation veloc- 
ities > 0.1 c. This finding strongly suggests that these 
sources are very young. Also a recent study of the inte- 
grated radio spectra of a number of CSS/GPS (Murgia et 
al. 1999) provides additional evidence for short lifetimes. 



It is nevertheless important to support these recent 
results with information on the ambient medium around 
these sources. 

The frustration scenario requires that an unusually 
dense ISM be present in CSS/GPS sources. This may con- 
tain a substantial cold phase and a large amount of dust, 
and, in order to stop the advancing jet, this medium has to 
have a large covering factor. The dust will therefore absorb 
and re-process a fraction of the optical and UV radiation 
from the AGN (and young stars, if present) higher than 
in the larger size, "non-frustrated" , radio sources, giving 
to the CSS/GPS sources an extra IR emission component 
in addition to the one from the disk/torus. 

Since this dense medium has to extend over several kpc 
from the nucleus (in order to stop also radio sources 10-15 
kpc in size), a large fraction of the dust will be moderately 
cold so that the emission is expected predominantly in the 
medium-far IR (MFIR) spectral region. 

By searching for emission from cold dust, we have 
therefore tried to obtain some constraints on this sce- 
nario. Since the emission peak is expected at relatively 
long wavelengths, ISO, with its capability of carrying out 
photometry at A > 60/im, (Kessler et al. 1996 ) was well 
suitable for such an investigation. 

This search is better suited to radio galaxies, rather 
than quasars, as in the former the contamination from 
the IR emission produced in the innermost regions by the 
AGN is lower by a factor 4-5 (Heckman et al. 1992). 



2 Throughout this paper we use Ho = 100 (Km/sec)/Mpc and 
qo=0.5. 



No definitive explanation, within the "unified models" 
scenario, has yet been given for such a difference (see 
also Sect. 6.5), but it is an observational effect that we 
have taken into consideration in planning the experiment. 
Moreover there are several ISO programs in the Guaran- 
teed Time allocation aimed at observing different samples 
of quasars, including a number of CSS/GPS, from which 
the relevant information could be retrieved. 

The project of which we present the results is the merg- 
ing of two independent programs (P.I.s C.F. and H.F.) 
which both got observing time in the first ISO Call for Pro- 
posals. A representative sample of CSS/GPS radio galax- 
ies, unbiased with respect to FIR emission, was selected 
from various samples of CSS/GPS radio sources. A sample 
of extended radio galaxies, matched in redshift and radio 
luminosity, was also selected for observations at the same 
wavelengths, in order to determine the level of MFIR emis- 
sion in extended radio galaxies as compared to CSS/GPS 
galaxies. 

The layout of the paper is the following. 

Sect. H describes the selection of the two samples. 

Sects. U and [| describe the ISO observations and the data 
reduction. In these sections we mention only what is rele- 
vant to the paper. For more details we refer the reader to 
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Table 1. CSS/GPS sample 













nominal observing time 


Source 


z 


LogP 2 .7GHz 


LS 


Sample 


60^tm 


90^m 


174/xm 


200/nn 






(W/Hz) 


kpc 




sec 


sec 


sec 


sec 


0108+38 


0.67 


26.79 


0.023 


PW 


256 


64 


128 


128 


3C49 


0.62 


26.89 


3.6 


3CR 


192 


64 


96 


128 


3C67 


0.31 


26.30 


6.8 


3CR 


256 


128 


— 


— 


0404+76 


0.60 


27.23 


0.53 


PW 


192 


96 


96 


192 


1244+49 


0.21 


25.63 


5.7 


San 


256 


64 


128 


128 


3C268.3 


0.37 


26.54 


3.9 


3CR 


256 


128 


32 


32 


1323+32 


0.37 


26.75 


0.18 


PW 


192 


64 


96 


96 


1358+62 


0.43 


26.81 


0.16 


PW 


64 


32 


32 


32 


3C303.1 


0.27 


25.60 


5.0 


3CR 


192 


96 


32 


32 


1607+26 


0.47 


26.92 


0.18 


PW 


64 


32 


32 


32 


1622+66 


0.20 


25.66° 


0.0013 


Sn 


64 


64 


64 


128 


3C343.1 


0.75 


27.24 


1.3 


3CR 






32 


32 


1819+39 


0.798 6 


27.23 


3.4 


PW 


192 


96 


96 


96 


1819+67 


0.22 


25.06 


0.02 


Sn 


192 


64 


32 


64 


1 Qon i on 
lozy+zy 


U.OU 


ZD.yo 


n i 
9.o 


T>WT 

L W 


64 


32 


64 


128 


2342+82 


0.74 


27.25 


0.66 


PW 


192 


64 


64 


128 


2352+49 


0.24 


26.18 


0.15 


PW 


64 


32 


64 


128 


median d 


0.43 


26.79 


1.0 






3C318 


1.574 c 


27.71 


8.2 


3CR 


64 


32 


128 


192 


1345+125 


0.12 


25.79 


0.11 


St 


64 


32 


32 


64 


1718-649 


0.014 


24.20 a 


0.0014 


Tin 


256 


64 


128 


128 


1934-63 


0.18 


26.65 


0.14 


Tzi 








64 


1946+70 


0.10 


25.00 


0.04 


Sn 


64 


64 


64 


128 



^3 GHz, and its P2.7GHZ has been extrapolated from the spectrum at high frequencies 



NOTES: 

a) the spectrum of this GPS peaks at 
6) revised redshift by Vermeulen et al. ( 1996) 

c) revised redshift by Willott et al. ( |2000| j 

d) the median values are for the abov e representative sa mple ( 0.2 < z < 0.8) 

Sampl e - PW: Pea cock & Wall ( |1982D ; 3CR: Laing et al. ( |1983 ); Sa n: San ghera ( [1995D ; Sn: Snellen et al. ( |l998| ); St: Stanghellini 
et al. Jl997|; |l998[); Tin: Tingay et al. (|l997[); Tzi: Tzioumis et al. (|l997|) 



the specific ISO and (IsoPhot Interactive Analysis) 

literature. 

Sects. |^ and H present the results and discuss them. 
Conclusions are summarized in bf Sect. ^ 

2. The sample 

Eighteen CSS/GPS radio galaxies in the redshift range 
0.2 < z < 0.8 and with radio size < 10 kpc were origi- 
nally selected. They are mostly from the 3CR catalogue 
at 178 MHz (Laing et al. 1983) and from the Peacock & 
Wall dl982| ) (P W) cat alogue of radio sources at 2.7 GHz, 
(see Fanti et al. 1995| ), but some have been selected also 
from other lists provided the selection was unbiased with 
respect to IR emission (see Table |l|). 

3 PIA is a joint development by the ESA Astrophysics Divi- 
sion and the ISOPHOT Consortium led by the Max Planck 
Institute for Astronomy (MPIA), Heidelberg. Contributing 
ISOPHOT Consortium institutes are DIAS, RAL, AIP, MPIK, 
and MPIA. 



Of these radio sources, 3C318 has recently been dis- 
covered to be a quasar at a redshift higher than originally 
believed (Willott et al. 2000| ) and therefore it has to be 
dropped from the sample, which counts then 17 CSS/GPS 
galaxies. However, for completeness, we have included in 
Table |B| also the data we obtained on this source. 

A sample of 16 3CR radio galaxies, with radio sizes > 
20 kpc and spanning a similar range in redshift and radio 
luminosity, was selected for comparison purposes. 

The two samples are presented in Tables & ||. 

Four additional GPSs, at redshifts < 0.2 were also se- 



lected for observations (see Sect. 5.3), although they are 
not going to be used in the discussion. They also are listed 
in Table [j], separate from the representative sample. 

3. Observations 

The sources were observed with ISOPHOT (Lemke et al. 



1996 ), sub-instruments C100 & C200, at the wavelengths 
of 60, 90, 174 & 200 /im. These wavelengths were chosen 
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Table 2. Control Sample 











nominal observing time 




Source 


z 


Log P2.7GHZ 


LS 


60/rni 


90^m 


174/xm 200Atm 






(W/Hz) 


kpc 


sec 


sec 


sec 


sec 


3C16 


0.41 


26.31 


135 


64 


32 


32 


64 


3C19 


0.48 


26.77 


20 








— 


64 


3C34 


0.69 


26.67 


174 








— 


64 


3C42 


0.40 


26.51 


85 








— 


64 


3C46 


0.44 


26.18 


533 


256 


64 


32 


64 


3C79 


0.26 


26.32 


212 


196 


32 


128 


128 


3C274.1 


0.42 


26.52 


470 


64 


32 


32 


64 


3G2( I .A 


U.I 1 


0^ no 
20.92 


Ol A 

214 


64 


32 


64 


128 


3C284 


0.24 


25.87 


409 


64 


32 


32 


32 


3C295 


0.46 


26.51 


20 








32 


3C299 


0.37 


26.44 


33 








32 


3C330 


0.55 


27.18 


192 








32 


3C337 


0.64 


26.96 


153 


64 


32 


32 


32 


3C401 


0.20 


26.12 


40 


196 


64 


32 


64 


3C441 


0.71 


27.03 


117 


64 


32 


32 


128 


3C459 


0.22 


26.18 


20 


64 


32 


64 


64 


median 


0.43 


26.51 


150 





in order to cover a spectral range as broad as possible 
with the best detectors. Actually some sources, mainly in 
the control sample, were observed in other programmes, 
including Guaranteed Time Observations, requiring us to 
observe them only at the wavelengths not planned in the 
other experiments (see Tables |l| & ||, where the wave- 
lengths missing in our programme are marked by a " — "). 

The CI 00 detector consists of a 3x3 matrix of pixels 
43.5" x 43.5" in size and C200 of a 2x2 matrix of 89.4" x 
89.4" pixels. Each pixel is an independent detector, which 
requires its own calibration. 

The diameter of the Airy disks of the Point Spread 
Functions (PSF) are given in Table ||. By comparison with 
the pixel sizes, it is clear that for C100 the Airy disk 
matches approximately one pixel size, while for C200 the 
PSF covers most of the matrix. Therefore detected sources 
will be visible mostly on the central pixel #5 with C100, 
while with C200 the four pixels should give approximately 
the same values (within the noise). Table || gives also the 
fraction, /psf, of light falling onto the central pixel (# 5) 
for C100 and on the four pixels for C200. 

Observations were made in chopper rectangular mode, 
with chopper throw 180" , in order to have, every few sec- 
onds (chopper plateau), an ON-source and an OFF-source 
measure of the same time length. The OFF data then have 
to be subtracted from the ON data, in order to extract the 
source signal. In this observing mode, the telescope points 
halfway between the source and the background positions, 
and a small mirror switches between the two sky positions. 
This introduces a vignetting error which is different for 
ON and OFF positions (see Sect. §). 

We applied for and obtained time also in the second 
"Call for Proposals" , mainly for the purpose of increas- 



Table 3. Airy disk size, e?Airy and /psf for the filters used 
in the observations 



filter 




/psf 


filter 


dAiry 


/psf 


/im 


["] 




fim 


["] 




C100 






C200 






60 


50.3 


0.69 


174 


134 


0.707" 


90 


75.5 


0.61 


200 


168 


0.794° 



a) Laureijs (1998) - private communication 



ing the observing time (originally planned on the basis 
of the ground estimates of the instrument performances). 
Therefore some of the sources were observed twice. 

The total ON-source time ranges from 32 to 256 sec, 
although the glitch removal and other options of the data 
reduction (Sect. 3) have shortened the nominal observing 
time. In Tables [l & || however, we give the total "nominal" 
ON-source time; for the sources observed twice we give the 
total. 



4. Data reduction 

The data analysis was carried on using mainly the PIA 
V7.2 software (Gabriel et al. 1997 ). Some further software 
has been kindly made available to us by Dr. M. Haas from 
MPIA (Heidelberg) or has been written by ourselves. 

PIA removes glitches due to cosmic rays, subtracts the 
dark current, corrects for drifts, derives the signals from 



Fanti et al. 



5 



the mmp^j, calibrates the data and makes corrections for 
vignetting. 

Special PIA features that we used are: 

1) We used the method of ramp subdivision, i.e. we di- 
vided each ramp into pieces 0.25 sec long (pseudo-ramps) 
to improve the accuracy of signal determination. 

2) to remove glitches we used the "two threshold glitch 
recognition" which is claimed to work better than the orig- 
inal one. In a number of cases, however, glitch "tails" re- 
main in the data. We wrote a simple IDL program to 
further statistically remove the glitch residuals. The noise 
is reduced by approximately a further 20%. 

3) Due to "memory effects" the chopper plateaus 
are not always flat, but in many cases show up- 
wards/downwards trends in the ON/OFF measurements. 
In absence of an appropriate program which fits these 
trends and extrapolates the asymptotic value, we used the 
quite crude PIA option which removes the first half of each 
chopper plateau, thus cutting in half the observing time. 

More detailed remarks on some other of these reduc- 
tion steps are discussed here below. 

4--1- Calibration 

The conversion Volt/sec into Watt is made by using an 
internal calibration source which provides the conversion 
factors for each of the matrix pixels (or detectors) . Then 
the conversion to Jy is made by using the known instru- 
ment and filter characteristics. The flux density scale may 
still be <J 20% uncertain at 60 and 90 /Ltm (see Sect. 5.2 ) 
and perhaps more at 174 and 200 fim. 

In addition, the matrices of the OFF measurements are 
not "flat" after PIA calibration, the pixel-to-pixel fluc- 
tuations being much larger than the instrumental noise. 
This is explained as a residual calibration error of the in- 
dividual pixels (flat fielding), that we estimate to be up to 
pa ±15%. This causes fluctuations from pixel to pixel of 
several tens of mJy, essentially due to the background, as 
the "true" signal from the sources is rather small. These 
errors, however, largely cancel out in the (ON - OFF) data 
and we have not attempted any correction. 

Finally, according to the ISO team, it seems that flux 
densities have to be increased by an amount which de- 
pends on the chopper frequency (i.e. the rate at which the 
chopper switches between ON and OFF positions) and on 
the signal difference between ON and OFF. At the time 
of writing this paper this correction has not been firmly 
established yet, although it is believed to be small for the 
C200 data. For the C100 data only the strong sources may 
be affected by this problem (see also Sect. p^2), while for 



4 The signal is read from the detector every 1/32 sec and ac- 
cumulated until the voltage reaches a saturation limit. At this 
time a "destructive" read-out occurs, which resets the integra- 
tion. The integral signal between two consecutive destructive 
read-outs is called a "ramp" and its slope gives the signal in 
Volt /sec. 



the average detections (see Sect. 5.4) we expect the cor- 
rection to be negligible. 

4-2. Vignetting corrections 

The vignetting corrections applied by PIA have an uncer- 
tainty of a few percent. As they are different in the ON 
and OFF measures, they do not cancel in computing (ON 
- OFF), as residual flat-field errors do. Given the typical 
values of the background (Table |j), this uncertainty intro- 
duces errors which are not negligible as compared to the 
source flux densities which are rather weak. 

At 60 and 90 /im, pixel # 5, where the source signal is 
largely concentrated, is not affected by vignetting and the 
source flux density is safely derived as explained in Sect. 
4.3 by dividing the signal on pixel # 5, by the /psf- 

At 174 and 200 /xm, where the source signal is spread 
over the four pixels, vignetting errors are a major cause of 
uncertainty, due to the higher background levels, not only 
as they increase the final noise, but also as they intro- 
duce systematic effects. Indeed, after application of PIA 
vignetting corrections, we got the peculiar result that the 
(ON - OFF) flux densities were, on average, systemati- 
cally negative, with a clear dependence on the background 
brightness. This effect is very likely due to residual vi- 
gnetting errors. 

We have therefore tried to estimate such corrections 
independently, from our own data. For this purpose we 
used a method suggested by Dr. M. Haas, from MPIA, 
whose application we take the responsibility for. We as- 
sumed that the pixels which are more central with respect 
to the telescope pointing, i.e. # 2 and #3 at ON, and # 1 
and #4 at OFF position, are, at the first order, unaffected 
by vignetting, while the other pixels are. 

As the PSF fills essentially all four pixels, they do see 
the same signal. Any deviation from this is attributed to 
vignetting. Considering only the sources with very strong 
background (^,20 MJy/ster), assuming that the source 
signal is negligible as compared to the background, and 
that the background does not change appreciably from 
ON to OFF, the four pixels at ON position ought to see 
approximately the same signal as the four pixels at OFF. 
Therefore the ratios of the signal in pixels # 1 and # 4 
ON (assumed to be the only ones affected by vignetting) 
to the signal of the corresponding pixels at OFF posi- 
tion (assumed to be unaffected by vignetting) give the vi- 
gnetting errors for these two pixels. The same procedure 
is applied to pixels # 2 and #3 at OFF position. Note 
that flat-fielding errors cancel out in taking the ratios of 
corresponding pixels. 

The correction factors deduced from the different 
sources are consistently in the range from 3 % to 5 % and 
were applied to all the sources. The estimated residual un- 
certainty is w 1.0 %. The application of these corrections 
eliminates the problem of the systematic negative (ON - 
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OFF) values and of their dependence on the background 
brightness. 

4-3. Flux density determination 

For C100 the source flux density should be derived by fit- 
ting a PSF of known width to the (ON - OFF) matrix 
values. In practice our signal is always very weak and will 
only be detectable on pixel #5. Therefore we have com- 
puted the source flux density simply by dividing the flux 
density falling onto pixel #5 by the appropriate value of 
/psf given in Table |[ 

For C200 the flux density has been obtained by sum- 
ming the (ON-OFF) four matrix pixels and dividing the 
result by the /psf in Table H. 

In the following we shall therefore use as units 
mjy /pixel to refer to the flux density falling onto one sin- 
gle pixel and mjy to indicate source flux density. 

4-4- Estimate of statistical flux density errors 

Besides the calibration and vignetting errors, additional 
uncertainty in the flux density measurements is caused by: 

a) instrumental noise (or simply noise, a n ), which includes 
all noise contributions along the signal chain (i.e. photon 
noise due to source and background, plus detector noise), 

b) confusion due to foreground galactic cirrus ( <t cc ) and c) 
confusion due to extragalactic background (c cgc ). Typical 
values are reported in Table |[ 

In the next paragraphs we discuss how they have been 
evaluated in the paper. 

a) Instrumental noise, <r n 

The nominal instrumental noise per pixel is given by 
the equation: 



NEP 2 



2NEP 2 ck + 2NEP? ec 



(1) 



where the NEP's are the Noise Equivalent Power of 
source, background and receiver, to be computed using 
parameters given in the ISO manuals and which depend 
on the ON-source time. The factors of 2 in Eq. (|l|) de- 
rive from the fact that the signal is obtained from (ON - 
OFF) measurements. ^From the values of the parameters 
necessary to compute the NEP's it is clear that, except 
for very high backgrounds (> 100 MJy/ster) or for very 
strong sources (j> 50-100 Jy), the contribution of the de- 
tector to the nominal noise is expected to dominate. 

Because the resulting performance of the instrument 
was worse than expected, we adopted a pragmatic method 
and estimated a' n from the data themselves, instead of 
using the parameters given in the manuals. 

To estimate a' n directly from the data, within each ob- 
servation we used the chopper sequence (i.e. the sequence 
of ON-OFF pairs of each source) to determine the statis- 
tics of the (ON - OFF) values. For C100 we used only pixel 



# 5, thus avoiding the residual flat-field (Sect. 4.1) and 
vignetting (Sect. 4.2) problems. For C200 we used the four 
pixels together, after application of vignetting corrections, 
and ignoring the residual flat-field corrections. 

The average (ON - OFF) gives the flux density within 
pixel # 5 for CI 00 and the flux density within any of the 
four pixels for C200. The r.m.s. distribution of the ON- 
OFF pairs gives the "actual instrumental noise" a' D (in 
mJy/pixcl) of that one observation. The software for this 
analysis is due to M.Haas (private communication). 

We find that a' n can change significantly from source to 
source, even for similar backgrounds and observing times. 
The median values, however, scale rather well with the 
inverse of the square root of the ON-source time. We do 
not find any difference as a function of the background 
brightness, which confirms that the detector noise is the 
dominant factor. The noise computed in this way is in 
agreement with that from Eq. [j] when adopting the pa- 
rameters of the ISO manual, except for A = 200/zm where 
our estimate is lower. This proves, to us, that no significant 
systematic error is present in our data and that glitches 
have been removed fairly well. 

Given the good behavior of the instrumental (mostly 
detector's) noise with l/Vt, we give in Table f|, for each 
observing wavelength, one single value for a' n , correspond- 
ing to an integration time of 64 sec. Note that with 
C200, since the four pixels are independent detectors, the 
flux density error due to the noise will be computed as 
cr n = 2 x a'Jf psf- 

b) Cirrus confusion noise, a cc 



This term has been evaluated by Gautier et al. ( 1992 ) 
who provide useful tables to compute er cc on any angu- 
lar scales (3 arcmin, in our case, the value of the chop- 
per throw) and also show that cirrus confusion depends 
strongly on the galactic background brightness, B, accord- 
ing to cr cc oc B 15±0 - 2 . 

Since our measurements include a contribution from 
the celestial fore/back-grounds and, at least at 60 and 90 
/im, a strong contribution from zodiacal light, in order to 
estimate, <r' cc (within a pixel), at all four wavelengths, for 
both C100 and C200 we used IRSKY software by IPAC 
which applies a model to remove the zodiacal light, pro- 
viding also the values for a pure galactic foreground. This 
estimate of o~' cc is accurate to within a factor of two. 

Note that, for A > 100/im, the foreground provided 
by IRSKY is extrapolated from the IRAS shorter wave- 
lengths, using the dust emission model by Desert et al. 
( 1990| ). This increases the uncertainty on a' cc . For consis- 
tency and to be conservative we generally used the IRAS 
background, which is, on average, higher than ISO's, ex- 
cept in those cases in which, by comparison with the other 
satellite values, the extrapolation from IRAS was obvi- 
ously wrong. In those cases we used the ISO background. 

We have also compared IRSKY values of a' cc with 
those reported in ISO manuals: while at 60 and 90 /zm 
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the two estimates do agree, at A = 174/xm the ISO man- 
ual estimate is over one order of magnitude greater than 
IRSKY's. We found no estimate at 200^m. 

As a reference, we report in Table |^ the values of a' cc at 
each wavelength, from IRSKY, for a foreground brightness 
typical of our observations. 

c) Extragalactic confusion noise, a cgc 

This is simply due to the piling up, within a beam, of 
faint galaxies not individually resolved by the instrument. 
The estimate of the magnitude of this effect, within each 
pixel, is reported in Table [I] and is derived by IRSKY, in 
the hypothesis of cosmological evolution of faint galaxies. 
These values are about a factor of two higher than in the 
case of a non-evolving model. 

Note that with C200, since the PSF covers most of 
the four pixels, we considered the cirrus and extragalactic 
confusion correlated across the matrix and therefore we 
adopted, as the flux density error, 4x (pixel error)//psF- 



Table 4. Noise, confusion and total flux density errors 



The total error on source flux density is then computed 



as: 



CTT 



cgc 



where a n is the individual source noise, a, 



cgc 



- . / /t — 4- ^ 4- „- (2) 

is the error 

due to the extragalactic confusion (which is the same in all 
fields at a given wavelength) and ct cc the cirrus confusion 
noise appropriate for the specific source background (zodi- 
acal light subtracted). In Table ^ total flux density errors, 
appropriate to the parameters in the Table, are also given 
as an example. 

As it can be seen from Table |], the extragalactic source 
and cirrus confusion (<J cgc and ct cc ) may be not negligible 
as compared to the instrumental noise and may limit the 
instrument sensitivity considerably or lead to spurious de- 
tections, as it will be further discussed in Sect. 5.1 . 

As an internal test, we used those sources observed 
twice (with different backgrounds) to check the goodness 
of error evaluation. 

At each wavelength, we compared the difference AS of 
the two measured flux densities with the combined total 
error (erf.) obtained by adding in quadrature the ctt of 
the two observations. If errors are estimated properly, the 
ratio AS/a^, computed for all the available pairs, should 
have zero average and standard deviation =1. 

In spite of the very poor statistics (8 sources observed 
twice with C100 and 4 with C200) the agreement is satis- 
factory. 

5. Results 

5.1. Individual Detections 

About half of the observed radio galaxies, both in the 
representative sample of CSS/GPS and in the comparison 
sample, have S > 3 x a n , at least at one wavelength. These 
could be considered "formal" detections, in the sense that 



A 


<r„(64 sec) 


/ 

°cc 


bckg 


/ 

^egc 


(Tt 


/jm 


mjy/pix 


mjy/pix 


MJy/sr 


mjy/pix 


mjy 


60 


27 


9 


(2) 


2 


42 


90 


12 


5 


(2) 


4 


22 


174 


23 


8 


(10) 


22 


148 


200 


61 


6 


(20) 


23 


195 



Values in parentheses are the background brightness at which 
a' cc is computed 



we got some signal out of the instrumental noise. This 
does not mean, however, that we have detected our target 
sources, since, as seen in Table |^, extragalactic and/or 
cirrus confusion could cause spurious detections: we might 
just have detected a cirrus filament or a fluctuation in the 
confusion. 

We have checked the reality of these "formal" detec- 
tions by looking for "negative" detections. We have only 
three such cases (see Tables [| and ||) suggesting that, af- 
ter all, in the "formal" detections some signal from the 
target is present. This is why we have kept them in Table 
||. The test for "claiming" a detection, however, is made 
by comparing S with ctt (Eq. |^) . 

We considered detected the sources for which S > 
5 x ctt at least at one wavelength (value in bold in Table ^j) 
and possibly detected those for which 3xctt<S'<5xctt 
(value in italic in Table |5|) . 

Table | gives the detected, possibly detected and 
"formally" detected sources, listing the flux densities at 
all four wavelengths (significant or not) and related er- 
rors. The sources detected at least at one wavelength are 
marked with **; those possibly detected are marked with 
*. In each column we give S, a n and ctt in mjy (1 a level). 

In Table |B| only two sources from the representative 
(1819+39) and control (3C459) samples are clearly de- 
tected at least at one wavelength. Nine more sources have 
flux densities which exceed 3 ctt, but are < 5 ctt, and are 
therefore possible detections. 

The quasar 3C318 is detected at a > 3ct level at both 
60 and 90 /im. 

We also detect clearly at all four wavelengths two of 
the four nearby GPS (1345+125 and 1718-64). 

In Table ^ we list the undetected sources with their S, 
CT n and ctt, in mjy. 

5.2. Comparison with IRAS data 

As stated in the IRAS Explanatory Supplement (1985) 
the sensitivity limits (3ct) of the all-sky survey for point 
sources (PSC) are quite high (~0.5 Jy at 60 fim and 1.5 
Jy at 100/xm) but can be sensibly reduced by co-adding 
s evera l scan s, as done, for instance by Heckman et al. 
( 1992j 1994 ), who obtained 3ct detections on individual 
objects of «100 and «400 mjy at 60 pLm and 100/xm re- 
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Table 5. Sources detected at different levels of reliability (mJy) 



Source 


A = 


60 /jm 




A = 


90 fita 




A = 


174 fim 


A = 


200 fim 




b 


On 




b 




CJ"x 


b 


On 


CJ"T 


b 


On 




* 


3C49 


42 


14 


15 


16 


17 


23 


335 


52 


146 


623 


102 


155 




1244+49 


22 


18 


21 


36 


16 


17 


187 


48 


133 


55 


134 


177 


* 


1323+32 


1 


16 


17 


45 


11 


13 


187 


36 


130 


306 


108 


157 




3C303.1 


-2 


18 


23 


45 


17 


19 


543 


117 


172 


368 


248 


274 




1622+663 


77 


30 


33 


43 


8 


12 


-85 


73 


145 


81 


89 


146 


** 


1819+39 


91 


19 


21 


147 


6 


12 


33 


53 


136 


8 


132 


176 




1829+29 


99 


40 


51 


129 


36 


80 


-741 


44 


313 


-348 


107 


171 


* 


2352+49 


8 


59 


61 


66 


27 


51 


204 


47 


343 


65,5 


130 


202 


* 


3C318 


118 


35 


36 


104 


21 


24 


277 


60 


150 


239 


68 


135 


** 


1345+125 


1470 


45 


46 


1171 


23 


24 


1273 


119 


172 


1471 


287 


310 


** 


1718-64 


499 


27 


33 


602 


13 


44 


2331 


41 


396 


2216 


77 


184 




1946+70 


56 


49 


51 


10 


11 


41 


-28 


52 


334 


433 


105 


184 




3C16 


-77 


33 


44 


-16 


25 


33 


181 


47 


164 


393 


156 


197 


* 


3C34 




















781 


149 


193 




3C46 


11 


17 


17 


32 


10 


13 


283 


115 


177 


212 


162 


200 


* 


3C79 


66 


24 


29 


37 


14 


37 


378 


29 


326 


625 


103 


181 


* 


3C277.2 


87 


25 


27 


-18 


13 


16 


289 


42 


133 


207 


152 


191 




3C284 


179 


47 


47 


18 


29 


30 


51 


150 


197 


86 


180 


214 




3C401 


2 


14 


17 


32 


18 


22 


238 


55 


158 


469 


213 


244 




3C441 


5 


28 


30 


-7 


19 


29 


11 


80 


246 


453 


120 


178 


** 


3C459 


698 


44 


51 


747 


45 


48 


668 


43 


157 


1164 


157 


196 



flux densities in bold characters are those for which S > 5 x or; 

flux densities in italic characters are those for which 3xo-t<S'<5xctt 

Table 6. Undetected Sources (mJy) 



Source 


A = 


60 


pirn 


A = 


90 


pirn 


A = 


174 




A = 


200 






S 


On 


Or 


S 


On 


or 


S 


O n 


Or 


S 


O n 


Ot 


0108+389 


9 


19 


19 


5 


11 


14 


45 


40 


149 


242 


106 


158 


3C67 


-2 


12 


18 


28 


13 


33 














0404+76 


7 


37 


39 


-14 


18 


33 


-645 


95 


286 


-217 


170 


220 


3C268.3 


17 


22 


24 


2 


15 


17 


34 


57 


137 


151 


224 


252 


1358+624 


54 


30 


31 


22 


29 


30 


124 


168 


210 


222 


213 


243 


1607+26 


71 


60 


62 


83 


43 


45 


-11 


122 


181 


302 


223 


252 


3C343.1 














-175 


84 


152 


-151 


102 


154 


1819+670 


23 


13 


16 


-26 


27 


29 


23 


86 


157 


-55 


199 


231 


2342+82 


-3 


30 


33 


32 


14 


36 


-68 


81 


329 





219 


265 


1934-63 








40 


93 


154 


3C19 




















630 


229 


256 


3C42 




















348 


132 


179 


3C274.1 


53 


49 


54 


-18 


32 


33 


6 


71 


144 


-237 


128 


173 


3C295 




















-111 


102 


154 


3C299 




















358 


148 


188 


3C330 




















191 


246 


272 


3C337 


-21 


33 


34 


43 


17 


18 


-141 


58 


138 


353 


192 


224 



spectively. The most sensitive IRAS observations were ob- 
tained by pointed obser vation s (see e.g. Neugebauer et al. 
1986 ; Golombek et al. 1988 ). The resulting noise varies 
from field to field and is typically in the range of 50- 
100 mJy (3cr) at 60^m and 100-300 mJy (3a) at 100/mi. 
These values have to be compared with our typical 3cr er- 



rors, given m Tables | and |, of «90 mJy at 60 /Lim and 
~75 mJy at 90/im (median values). Therefore, while at 
A = 60^m ISO is not much more sensitive than IRAS, at 
A = 90/im the situation is improved. 

We searched the literature to see whether any of the 
sources in our sample has been detected by IRAS at any 
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Table 7. Available IRAS flux densities (mJy) 



Source 


60/im 


WOum 


ref. 


3C79 


179±9 


<330 


1 


1345+125 


2098 


1738 


2 


3C318 


170±26 


300±70 


1 


3C459 


683 


708 


2 



ref: 1) Heckman et al. (1992); 2) Golombek et al. (1988) 



wavelength. The galaxies from our lists detected by IRAS 
at 60 and/or 100 /im are reported in Table |[ 

On the strongest source (1345+125) about 30% of the 
IRAS flux density is missing in our observations. This is 
the case where chopper frequency corrections, that we 



have not applied (Sect. 4.3), are expected to be larger. 
For 3C459, which is not as strong, the agreement with the 
IRAS data is rather good. We are therefore confident that 
the missing correction is definitely minor for the fainter 
sources and negligible for the median flux densities that 
we report. For the fainter sources, we note that 3C79 is 
possibly detected in our observations only at wavelengths 
longer than 100/im and that 3C318 at 90 jiva. is fainter 
than in IRAS, although the two flux densities are still 
within ~ 2.5 of the combined errors. 

Conversely 3C284 possibly detected (3.8 a level) at 
60 /im in our ISO observations, has only an upper limit in 
IRAS (Impey & Gregorini 1993| ) not too different from our 
measurement. The remaining sources detected (1819+39) 
or possibly detected at 60 and/or 90 ^m by ISO, are 
weaker than the IRAS upper limits (Heckman et al. 1994[ ). 



5.3. Comments on some individual sources 
1345+125. 

This radio source is identified with a galaxy with "two 
nuclei" , of which the western one has an optical spectrum 
consistent with Narrow Line radio galaxies and Seyfert 
type II (Gilmore & Shaw 1986). There is some confu- 
sion in the literature about which of the two hosts the 



radio source. Gilmore & Shaw (1986) and Baum et al. 
( |1988| ) indicated the eastern (non Seyfert nucleus). More 
recently, however, Capetti (1998, private communication), 
Stanghellini et al. (1998) and Evans et al. ( 1999| ) presented 
convincing evidence for the association of the GPS source 
with the western Seyfert nucleus. This double nucleus 
galaxy is considered a merger which has triggered the ra- 



dio activity (e.g. Heckman et al. 1986; Sanders et al. 1988; 
Surace et al. 1998] ). Although the optical spectrum is char- 
acterized by Narrow Lines, broad Paa (FWHM ~ 2600 



km s 1 ) is seen in th e nea r-IR (Veilleux et al. 1997 ). In 
addition, Hurt et al. ( |l999|) report HST/FOC ultraviolet 



observations which show extended polarized UV contin- 
uum light. Both the UV polarization and the broad Paa 
indicate that 1345+125 has a hidden (but not too deeply 
hidden) quasar. This is the first GPS radio galaxy de- 
tected in X-rays (£x(2-iokcv) ~ 2 x 10 43 erg s -1 , O'Dea 
et al. |200C| ) . 

We estimate a temperature of ~ 30 K and a dust mass 
of ~ 3.5 x 10 s M Q for the component which makes the 
dominant contribution to our ISO data (Table [To]). Al- 
though our ISO observations are not sensitive to gas at 
temperatures below about 25 K, the CO observations are 
sensitive to gas with kinetic temperatures to within a fac- 
tor of a few of the microwave background temperature 
(e.g. O'Dea et al. |1994| ; Allen et al. |1995[ ). For a stan- 
dard Galactic gas/dust ratio of 100 this corresponds to a 
mass of molecular gas ~ 3.5 x 10 10 M Q which is in good 
agreement with the estimate derived from CO observa- 
10 10 M n , Mirabel et al. Il989l Evans et al. 



tions (~ 3 x 10 iu M , Mirabel et al. _ 
1999D - This agreement between the ISO and CO estimates 
of the molecular gas mass may be fortuitous since the gas- 
to-dust ratio may vary between galaxies and in general is 



not well determined (e.g. Goudfrooij et al. 1994; Falco et 
al. p999h. 



3C 318 

This object was classified as an N-galaxy at a redshift 0.75 
by Spinrad & Smith (197£) and considered to be a Nar- 
row Line (NL) radio galaxy by Hes et al. 1995. Recently 
Willott et al. ( 2000| ), from UKIRT spectroscopic observa- 
tions, detected broad H Q and revised the redshift at 1.574. 
The object is now classified as a quasar. It is a strong 
X-ray emitter (Taylor et al. 



dated for the new redshift, 

c .-i 



of L 



1992) with a luminosity, up 
8 x 10 44 



X(0.5-4.5Kcv 



erg 



1622+663: 

This object has a broad a H Q line. It can be classified as 
a Seyfert 1.5 (Snellen et al. 1999). Its broad-line region 



seems to be heavily obscured. It is marginally detected at 
90 ^m. 

1718-649 

The source is identified with the 12.6 blue magnitude 



galaxy NGC 6328 (Savage 1976), whose nucleus has a 
LINER type optical spectrum (Fosbury et al. |l977p . The 
galaxy is quite peculiar, having the appearance of a high 
luminosity elliptical with faint outer spiral structure and 
containing a large amount of atomic hydrogen (Veron- 
Cetty et al. 



1995) 



The object is regarded as a merger 
of two galaxies, at least one of which is a gas-rich spi- 
ral, in the process of forming an elliptical. A very strong 
source (a star?) is listed in the IRAS PSC, at 3.7 arcmin 
from 1718-649. Inspection of IRAS maps at 60 and 100 
fim shows a fainter extension of this bright source to the 
North, at the position of 1718-649. 
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FIR properties of CSS/GPS radio galaxies 



The dominant component of dust contributing to our 
ISO detection has an estimated temperature of 20 K and 
a mass of ~ 0.1 x 10 s M (Table |o|) which would cor- 
respond to a gas mass of ~ 10 9 M Q . This is at least an 
order of magnitude less than the mass of HI (~ 3.1 x 10 10 
M Q , Veron-Cetty et al. |1995[ ). The HI is distributed over 
the entire galaxy and much of it lies in a partial ring with 
a diameter of 95" (19 kpc). Possibly, much of the dust in 
this object lies at large distances from the nucleus and is 
relatively cold, thus escaping detection in our ISO obser- 
vations. CO observations would be useful. 

1819+39 

Its redshift has been re-measured by Vermeulen et al. 
( |1996| ) who also show that this object has broad H^. Thus, 
this object appears to be Broad Line radio galaxy. 

1934-638 

The radio source is associated with the brighter galaxy of 
a pair of compact galaxies (ss 8 kp c sepa ration) sharing 
a com mon envelope (Jauncey et al. 1986 ). Fosbury et al. 
( |1987| ) suggest larger amounts of gas and dust than gen- 
erally found in typical Narrow Line radio galaxies. Tad- 
hunter et al. ( 1994 ) found presence of polarized light, sug- 
gesting an anomalous environment (see also Morganti et 
al. |1997| ). It is undetected at the only observed wavelength 
(200 urn). 

1946+70 

This is one of the th ree GPS sources from the sample of 
Snellen et al. ( 1998|) , and it is much weaker than those 
from the 3CR and PW samples. It is "formally" detected 
at 200/im. 

3C459 

The optical object has been classified as N-galaxy by Spin- 
rad et al. (1985). It was studied in detail at radio and op- 



Table 8. Average flux densities for the CSS/GPS and the 
comparison sample 



tical frequencies by Ulvestad ( 1985Q , who concluded that 
this galaxy may have undergone a lar ge amount of fairly 
recent star formation. Hes et al. ( 1995 ) classified it as Nar- 
row Line radio galaxy and suggest that a fraction of the 
FIR emission may be due to heating of the dust by the 
moderately young blue stars. 

5.4- Average Detections 

We have looked for "statistical" emission from both the 
representative CSS/GPS sample (i.e. for 0.2 < z < 0.8) 
and the control sample. We have computed both the mean 
and the median flux density, although in the following we 
use the second one, which is considered to be more robust 
since it is less affected by extreme flux density values. In 
computing the mean flux density we excluded 3C459 for 
C100, since its large detected flux density would affect 
the results too much (the source was used, instead, in 
computing the median). 



A 
fim 


CSS/GPS 
mean median 
mjy 


Control 
mean median 
mjy 


60 


36 


± 


24 


201- 


34 


± 


48 


oz -30 


90 


39 


± 


22 




11 


± 


18 


2 r + 15 


174 


-2 


± 


156 


04+90 


198 


± 


146 


2 l Q+90 
zlu -195 


200 


140 


± 


130 


iio -116 


370 


± 


168 


qcc + 120 



errors on both mean and median are 2 a (see text) 



Values are given in Table |8J. The uncertainty on the 
mean and median are derived from the flux density dis- 
tribution itself. Note that, as the error distribution of the 
median is not Gaussian, we choose to give the 95% confi- 
dence uncertainty, equivalent to a 2 a limit for a Gaussian 
distribution. To make the comparison easier we then give 
2<T errors also on the mean. 

At 60 and 90 /im there is clearly a 3 ex detection for 
the CSS/GPS sample. In the comparison sample, where 
the statistics is poorer, as we do not yet have data for 
w 1/3 of the objects (see Table ||) , both mean and median 
flux densities, compared with the errors, do not indicate 
a significant detection. They are, however, largely consis- 
tent with the CSS/GPS results. This is re i nforc ed by a 
comparison done with Hcckman et al. ( 1992 ; 1994 ). These 
authors give co-added flux densities for a sample of radio 
galaxies. In their "near sub-sample" (0.3 < z <^ 0.85, me- 
dian z — 0.48), composed of 5 CSS/GPS and 36 large size 
sources, and therefore similar to our comparison sample, 
the median flux densities (Sen ~ 25+5 mjy, S100 ~ 40+12) 
are consistent with ours, but have smaller errors, due to 
the much larger number of sources in their sample. 

At 174^m and 200/im radio galaxies in the control 
sample seem stronger FIR emitters, although the differ- 
ence is not significant at either wavelength. Averaging the 
two wavelengths the difference is around 2a. 

We have also applied a Kolmogoroff-Smirnov test, 
which shows that the hypothesis that the two galaxy sam- 
ples are extracted from the same population is acceptable 
at a 99% confidence level at all wavelengths. 

Considering also the results from Heckman et al. ( |1992| ; 
1994), we estimate that any IR excess 



for CSS/GPS 

sources does not exceed (at a 95% confidence level) ~ 25 
mjy in the range 60-90 /im and ^100 mjy in the range 
170-200 /im. 

Our first conclusion then is that we do not have any 
statistical evidence that the MFIR properties of the galax- 
ies associated with CSS/GPS sources are much different 
from those of the galaxies associated with large size radio 
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Table 9. Average flux densities and FIR luminosities for 
the whole sample 



A 


mean 

mjy 


median 


Lfir/L 


Lfir 
erg s _1 


60 
90 


35 ± 22 
29 ± 16 


20^ 
09+15 

oz -16 


6 x 10 10 


« 2.4 x 10 44 


174 
200 


74 ± 114 
256 ± 112 


48ijf 
909+123 

ZJZ -114 


« 10 11 


r;5x 10 44 



errors on both mean and median are 2 a 



sources (see also Sect. 6.2 ) consistent with the results of 
Hcckman et al. ( 1994 ). The small and large sources we 
observed seem to have similar average luminosities at all 
four wavelengths, and hence a similar spectrum between 
60 and 200 /im (observer frame). This will allow us to 
combine the two sets of data in the discussion (Sect. |6j~ 

6. Discussion 

6.1. FIR tail of synchrotron emission 

We have collected all the flux densities available in the lit- 
erature at radio wavelengths for every source in the sam- 
ple. No evidence is seen for bright flat spectrum cores at 
the highest available frequencies (often up to 230 GHz, 
1300 /im). The spectra can be well fitted either with a 
power law or with a model which includes spectral cur- 



vature due to synchrotron losses (see Murgia et al. 1999, 
for details). For each source, using the fitted spectrum, we 
have extrapolated the flux densities at 60, 90, 174 and 200 
yum. The extrapolated average synchrotron flux density is 
very small (<C, 2 mjy at 200 /im and less at the other wave- 
lengths) and does not contribute significantly to the FIR 
emission. Therefore we attribute the emission we find to 
genuine thermal radiation. 

6. 2. FIR luminosities of the combined CSS/ GPS and 
comparison sample 



As discussed in Sect. 5.4, our data do not show any sta- 
tistically significant evidence for differences in the FIR lu- 
minosities of CSS/GPS and large size radio galaxies at 60 
and 90 /im, the situation being somewhat more uncertain 
in the 174-200 /on interval. 

Now we combine the two samples together into a single 
sample of powerful radio galaxies at intermediate redshift 
to be studied on its own, independent of the type of asso- 
ciated radio source. In Table ^, we report the mean and 
median FIR flux densities of the combined sample and the 
average FIR luminosities in the 60-90/im and 174-200 /im 
ranges. 
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Fig. 1. Radio versus IR luminosity. Shaded areas repre- 
sent the regions occupied by spirals (below) and low power 
radio galaxies (above) from Knapp et al. ( 1990 ). The tri- 
angle marks the median values for the present sample of 
radio galaxies, while the empty and crossed squares are 
the low— z and high— z quasars respectively (Heckman et 
al. [I992D 



The FIR luminosities from the median values of the 
whole sample are calculated as: 

where S v is the median flux density at the frequency v 
and -Dl is the luminosity distance. 

The average Lfir are in the range (0.6 — 1.0) x 10 11 
L Q , or (2-5) X 10 44 erg s -1 . The brightest objects, 1819+39 
and 3C459 , which are several times more luminous than 
average (Table [l(]), have peculiarities, as commented in 
Sect. |5.3| , which may well account for their larger FIR 
luminosities. 

These luminosities are much higher than those found 
in nearby radio quiet elliptical galaxies, which are in the 
range of up to m 1 O 8 L (Goudfrooij & de Jong 1995 
Bregman et al 



1998) 



We point out that the median FIR luminosity we find 
follows the extrapolation of the Lfir vs L ra dio plot for 
nearby radio galaxies presented in Knapp et al. (1990) 
(see Fig. |). 

6.3. Temperature and Dust masses 

On the hypothesis that the dust is transparent at these 
wavelengths, the MFIR radiation due to dust at a uniform 
temperature T follows a modified Planck law given by 



S(u Q ) 



(x(i/ e )B h . h .(T, v e )M dust (l + z) 



Dl 



(3) 



where vq and v e are the observed and emitted frequencies, 
/i is the absorption coefficient, -Bb.b. is the brightness of 
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Table 10. Luminosities and dust masses for the bright sources 



Source 


Log(L 
60-90/mi 


fir/L©) 
174-200/im 


Hot Component 
T(K) Mass(M ) 


Cold Component 
T(K) Mass(M Q ) 


1345+125 


11.6 


12.2 


70 


4.4 xl0 b 


30 


3.5x10 s 


3C 318 


12.1 




100 


8 xlO 6 






1718-64 


9.3 


9.5 


60 


2.9xl0 4 


20 


0.1 xlO 8 


1819+39 


12.2 




65 


1.5xl0 7 






3C 459 


11.9 


11.5 


70 


4.6 xlO 6 


25 


3.7 xlO 8 



1000 



36 000 — 2000 GHz), can be described by a power law: 
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Fig. 2. IR spectrum of radio galaxies (crosses: IRAS, di- 
amonds: ISO) and quasars from Heckman et al. (1992) 
(asterisks). Error bars are 1 a 



a black body at temperature T, .Mdust is the dust mass, 
Z?l is the luminosity distance and z is the source redshift. 
The absorption coefficient is rather uncertain. In order 
to keep consistency with the values generally assumed in 
the literature (e.g.: Hildebrand 1983) ; Knapp et al. 1990| ), 
we take fi(X) = 10 x (250/A(/im)) for A < 250/xm. With 
respect to other assumed values, this tends to overestimate 
the dust mass (e.g. by a factor of two as compared to the 
n(X) assumed by Hughes et al. 



1997) 



The MFIR average spectrum for our sample is shown 
in Fig. |[ In order to describe it better, we used also the 
following additional upper limits or median flux densities: 
< 5 mjy at 1.3 mm (230 GHz), as maximum flux density 
excess over the extrapolated synchrotron spectrum (Mur- 
gia et al. |1999|) ; < 9 mjy a t 12/im and (12 ± 4) mjy at 25 
/jm (Heckman et al. 1992 near sub-sample). 

The dependence on frequency of the FIR flux den- 
sity in the observed wavelength range 12-200 /im (f re st ~ 



At lower frequencies, however, the spectrum of the FIR 
component must turn over to meet the 230 GHz upper 
limit. 

This spectrum is clearly not consistent with a single 
temperature. A multi-temperature fit (as in models by 
Sanders et al. 198£) would give a plausible range of tem- 
peratures between ss 25 K and > 120 K in the wave- 
length range we have considered here. However, as the 
data quality is not sufficient to allow a meaningful multi- 
temperature fit, we have taken a simple two tempera- 
ture approximation, with Ti = 80 K and T 2 = 25 K. 
With corresponding dust masses M\ w 5 x 10 5 M Q and 
M 2 « 2x 10 8 M Q , we can reproduce the FIR flux densities 
in the range from 60 ^m to 200 /im. 

Any dust component colder than 25 K is poorly con- 
strained by our data. For instance, up to 10 9 Mq would be 
consistent with the 174/200 /im median flux densities and 
the 1.3 mm upper limit if the temperature were < 20 K. In 
order to probe such low temperatures, measures at longer 
wavel ength s would be required, as with SCUBA (Holland 
et al., 1999 ). To our knowledge, such data do not exist for 
our sources. 

A comparison of our average mass estimates with those 



for nearby radio quiet ellipticals (Bregman etal. 1998) 
and nearby radio galaxies (Knapp et al. 199C ) is not 



straightforward, as these estimates depend on the wave- 
lengths avai lable and on the assumed temperatures. Breg- 
man et al. ( 1998 ) evaluate the temperatures from the 60 
/xm to 100/im flux density ratio (« 30 K) an d obtain 
10 5 < Mdust / Mq < 10 7 . Knapp et al. (|199C|) assume 



T w 19 K and derive 10 6 < M dust /M Q < 10 8 . 

In our opinion, it is far from being excluded that the 
range of dust masses in these samples and in ours are 
similar, and the different estimates may be due to either 
the different assumed temperatures and/or to the different 
wavelength ranges explored. 
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Using the two temperature approximation, we have 
also estimated dust masses for the brighter sources de- 
tected, or possibly detected, at more than one wavelength. 

The results of this calculations are presented in Table 
|To| . Compared to the median values of the total sample, 
the sources 1345+145, 3C318, 1819+39 and 3C459 seem 
to have a higher mass of the warm dust component, while 
the cold dust masses seem to be similar. As the first three 
objects are a Seyfert type galaxy, a quasar and a Broad 
Line radio galaxy, we may expect that a favorable orienta- 
tion of the objects with respect to the line of sight allows 
us to see the warmer dust associated with the accretion 
disk. 



6.4- What heats the dust? 

According to the unification models, the radio galaxies in 
our sample, which are very powerful in the radio band 
and have large luminosities in the Narrow Emission Lines 
(Loin « 10 43 erg s-\ Gelderman & Whittle, |1994| ), are ex- 
pected to harbor an obscured quasar (see, e.g., Falcke et al. 
1995| , in the context of " jet /disk-symbiosis"). A large frac- 



tion of the MFIR emission would then originate from the 
dust in the absorbing torus itself, heated by the UV/X ra- 
diation emitted at its center, with some contribution from 
the interstellar medium, still heated by the central AGN 
but now distributed over a scale of several kpc, and there- 
fore with temperatures decreasing outwards. We expect 
this interstellar dust to emit its radiation at the longer 
wavelengths. 

If an obscured AGN is present, adopting a heating 
model as in Sanders et al. ( |1989| ), we expect dust tem- 
peratures in a range similar to that estimated in Sect. 6.3. 

If, instead, no AGN is present and heating is due to 
the stellar radiation field only, the dust temperature could 
be <^ 30 K, as in nearby radio quiet elliptical galaxies 
(Bregman et al. 1998| ). The emission would then peak at 
A > 200/im. An additional contribution to heating could 
arise from a burst of star formation, as may be the case 
in 3C459. 

The interpretation based on heating by an obscured 
quasar is supported by the correlation between radio emis- 
sion and MFIR emission (Knapp et al. 1990; Heckman et 
al. |1994 ), with which our data are also consistent (Fig. 
[l]). Also the strong luminosity in Narrow Lines, which is 
often interpreted as powered by photo-ionization, gives 
support to this interpretation. In addition we note that 
the MFIR average spectrum (for A rcs t > 20/mi) in nearby 
radio quiet ellipticals (Bregman et al. |1998| ) is definitely 
steeper than that of our radio galaxies and of those in 
Knapp et al. ( |1990p . The spectral index of the former is 
chrqe — 2.3 ± 0.2, as compared to a us ~ 1.0 ± 0.2 in our 
sampl e (Sect. |6^2| ) and aKnapp — 1-2 ± 0.2 in Knapp et al. 
( |l990| ). This suggests that the heating processes are dif- 
ferent, with the one in radio galaxies being more effective. 



We conclude that our observations are consistent with 
the MFIR luminosity in these radio galaxies being pow- 
ered primarily by the AGN, and assume that most of the 
radiation comes from the "obscuring torus" . 

6. 5. Implications for the circumnuclear torus 

Our average FIR luminosities, consistent with the results 
of Heckman et al. (1992), are about a factor 4-5 lower 



than the FIR luminosities of radio quasars of comparable 
radio luminosity (see Fig. |2|). 

This causes problems for the unification models (see 
Heckman et al. 1992), since, if the torus dust is transpar- 
ent at FIR wavelengths, its emission should be orientation 
independent and the quasars should be as luminous in the 
FIR as the radio galaxies. On the other hand, if the disk is 
optically thick in some wavelength range, one should jus- 
tify the constant flux density ratio between radio galaxies 
and quasars over a broad wavelength range (from 12 to 
100 Aim). 



Hes et al. (1995) discuss the possibility that a beamed 
non-thermal component is present in quasars and con- 
tributes at the FIR wavelengths. 

6. 6. Back to the frustration scenario for CSS/ GPS 
sources 

The original aim of this work was to study the cold phase 
of the interstellar medium in radio galaxies, in order to see 
whether its density is high enough to support the "frus- 
tration scenario" for CSS/GPS sources. In this scenario, 
CSS/GPS and large size radio galaxies have similar life- 
times, so that the difference in size (more than a factor of 
100, on average, in our two samples, Tables [J and ^) re- 
quires an environment significantly denser for the former 



(De Young 1993; Fanti et al. 1995) 



The hot component of the ISM is too tenuous and has 
insufficient pressure to confine the radio source (e.g. O'Dea 

1996D - Thus the radio source 



et al. 1996; Readhead et al. 



must be confined via interaction with relatively cold dense 
material along the path of propagation. However, as antic- 
ipated in the introduction, even if a lot of dense cold gas is 
present in the host galaxy, distributed in the disk/torus of 
the "unified models" (which is perpendicular to the radio 
axis) it will have no effect on the radio source propagation. 
"Frustration", if present, must be due to the diffuse cold 
phase of the interstellar gas. 

/From ram pressure computations, the estimates of the 
density of the interstellar gas required to frustrate CSS ra- 
dio galaxies imply masses > 2 x 10 10 M within a volume 
of « 10 kpc in radius (Fanti et al. 1995), corresponding, 
for a gas to dust ratio as in our own Galaxy, to an excess 
dust mass > 2 x 10 8 M Q . This would produce an excess 
in the FIR emission of the CSS/GPS, as compared to the 
large size radio galaxies, the exact amount depending on 
the dust temperature. The upper limits to any IR excess 
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derived in Sect. 5.4 convert into an excess dust mass of 



<; 0.6 10 for T ;> 30 K, (if heating is due to a power- 
ful AGN), which is not enough for the frustration model. 
Even if we attribute the whole MFIR emission we detect 
from CSS/GPS to a frustrating medium only, and not to 
the dusty disk/torus, the implied mass (for T J£ 30 K), 
would be still not enough, the upper limits to the differ- 
ences in the MFIR emission being similar to the actual 
flux densities of CSS/GPS sources. 

As a consequence, we may be confident that there is no 
evidence for a component of the interstellar medium, with 
temperatures ^,30 K, significantly denser or more massive 
in CSS/GPS galaxies to cause frustration. 

For lower temperatures, our FIR observations are not 
sensitive enough to firmly constrain larger masses, since 
the emission would peak at A > 200 /im, and the uncer- 
tainties we have on the flux density limits at long wave- 
lengths are too large to provide safe constraints. For in- 
stance, if the temperature were «20 K, up to 5xlO 8 M 
of dust would be permitted. This should be probed by ob- 
servations at longer IR wavelengths (e.g. with SCUBA). 
An additional strong constraint on such large dust masses 
is obtained from the color properties of the host galax- 
ies. Dust masses as required by frustration, spread over ~ 
10 kpc, would produce absorption in the optical band. 
The implied hydrogen column density, iVn > 10 22 g 
cm~ 2 , corresponds to an optical depth in the visual r v « 
7Vh/(2 10 21 ). Applying a model of homogeneously mixed 
dust and stars (as in Goudfrooij & de Jong 1995 ), the 
above dust mass implies A v « 2 mag., and color redden- 
ing E(B-V) ~ 0.7, which is inconsistent with the data of 
de Vries et al. fll998|). 



If "frustration" occurs only for GPS (sizes < 1 kpc), 
the implied masses are lower, but their temperature woul d 
be larger (> 50 K, for a model as in Sanders et al. 1989| ), 
due to the closer proximity to the hidden quasar, so that 
the expected IR luminosity would be even more discrepant 
with what we observe. 

The present conclusion is therefore that the FIR radi- 
ation we detect with ISO does not imply an amount of gas 
(heated by the AGN) large enough to support the frustra- 
tion scenario of CSS/GPS radio sources. This is consistent 
with the results of Owsianik et al. (1998) and Owsianik 
& Conway (1998), who found fast proper motions of the 
outer lobes in some small symmetric GPS, and with the ra- 
dio spectra analysis of CSS/GPS by Murgia et al. ( |1999| ). 



We finally note that recent observations with SCUBA 



(Archibald et al. 2000) of high redshift radio galaxies, both 
CSS and large size ones, do not show any significant dif- 
ference between the two classes of sources. 

The MFIR properties of CSS/GPS, however, point out 
a new problem for the alternative interpretation. An im- 
portant ingredient of the "youth scenario" is that the 
radio luminosity is required to decrease with increasing 



source size, in order to explain the numerical proportion 
of CSS/GPS to large size sources. In the current models, 
this luminosity evolution is essentially due to expansion 
of the lobe volumes in an external medium of decreasing 
density and not to a decline of the power carried by the jet 
with time. CSS/GPS sources would evolve, with increas- 
ing size, toward lower luminosity larger size radio sources. 
This evolution in luminosity may also be described in 
terms of a higher efficiency for converting the beam power 
into r adio band radiation (see also Gopal Krishna & Wiita 



1991 ) in the early phases of the source evolution. 



According to this idea, our two samples (CSS/GPS 
and comparison) , which have very different sizes but simi- 
lar radio luminosities, would have "radio engines" of power 
different by a factor 10 or so (CSS/GPS being less pow- 
erful). As in the past there have been strong suggestions 
that the jet power is proportional to the bolometric lumi- 



nosity of the AGN (e.g. Baum & Heckman 1989; Rawlings 
& Saunders |l99l| ; Falcke et al. |1995| ), one would expect the 
CSS/GPS to be definitely less powerful in the MFIR (and 
in the NL emission). This is not seen in our observations, 
at least at 60 and 90 /im. 

Although we have no obvious solution to this prob- 
lem, we speculate here about a couple of possibilities. 
First there may be an additional source of heating arising 
from conversion of a fraction of the jet power where the 
jet interacts with the interstellar medium. This may be 
more effective at distances close to the nucleus and less 
so further out as the source ages. The existence of emis- 
sion line gas aligned with the radio source in low redshift 
CSS sources is consistent with strong interaction between 
GPS/CSS source s and the dense clouds in the ISM (de 
Vries et al. 1999 ). Such an effect should roughly compen- 
sate for the lower radiation power from the central engine 
in CSS/GPS. 

The second possibility is related to the scenario where 
the radio emission is triggered by a merger event. We spec- 
ulate that in the early phase of the radio source life a large 
fraction of gas/dust is not fully settled in a disk/torus 
but is still distributed over a much larger solid angle, as 
seen from the central continuum source. In this scenario 
the fraction of intercepted and reprocessed UV radiation 
could be much higher in the younger/smaller sources, in- 
creasing their MFIR emission. We might also expect the 
additional obscuration in the smaller sources to result in 
lower emission line luminosities and/or redder colors. This 
is not generally seen in GPS and CSS sources, though 
there are a few ca ses which are c onsis tent with this pic- 
ture (O'Dea |i"998[ de Vries et al. |l998|) . In addition, this 
hypothesis would require the gas/dust to settle into a disk 
on a timescale much less than the age of a radio source 
(10 7 ~ 8 yr) in order to remove the extra obscuration before 
the sources propagate to scales larger than tens of kpc. 

At present, both of these explanations appear rather 
"ad hoc" . Thus, this issue remains a problem for the cur- 
rent evolution models. 
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7. Conclusions 

a) We have presented ISOPHOT observations at A = 
60, 90, 174 and 200 of CSS/GPS radio galaxies and of 
a matched comparison sample of extended radio galaxies. 
A minority of objects are detected individually at one or 
more wavelengths in both samples. 

b) We have co-added the data for each sample to ob- 
tain mean and median flux densities. The extrapolated 
radio spectrum under-predicts the observed MFIR flux 
densities, arguing that the MFIR is due to dust rather 
than to synchrotron emission. 

c) We find no significant differences in the MFIR flux 
densities of the two samples. Our results are then consis- 
tent with the CSS/GPS and the extended radio galaxies 
having similar MFIR luminosities. For dust temperatures 
^,30 K, the deduced masses of the interstellar gas are lower 
than required by the frustration scenario. For lower tem- 
peratures larger masses would be allowed for by our data, 
but they would produce a large amount of obscuration and 
reddening in the optical, which is not seen in the existing 
data. 

All this argues against the CSS/GPS sources being 
"frustrated" by a dense ambient medium. 

d) Since no significant difference is seen between the 
two samples, we have combined them for further analysis. 
The average Lfir is in the range (0.6 — 1.0) x 10 11 Lq or 
(2 — 5) x 10 44 erg s _1 . Over the wavelength range of our 
observations, the spectrum can be fitted by a single power 
law with a ~ 1.0 ± 0.2. 

e) We have fitted simple two temperature (T\ = 80 
K and T 2 = 25 K) models to the IR spectrum and 
have derived dust masses of Mi S3 5 x 10 5 M Q and 
M 2 « 2 x 1O 8 M . The mass of the cold dust appears higher 
than found in radio quiet elliptical galaxies, although this 
difference may be due to the differences in the assumed 
temperatures and the sampled rest-frame wavelength cov- 
erage. 

f) Our observations are consistent with the MFIR lu- 
minosity in these powerful radio galaxies being mostly 
powered by an obscured AGN, although in some objects 
a contribution from star formation is possible. 
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